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Abstract 


An approximate theoretical method is 
presented which determines the limit cycle 
behavior of a helicopter model which has 
one or two nonlinear dampers. The 
relationship during unstable ground 
resonance oscillations between lagging 
motion of the blades and fuselage motion 
is discussed. An experiment has been 
carried out on using a helicopter scale 
model. The experimental results agree with 
those of the theoretical analysis. 
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geometrical length of the mode, 
see Fig.l 

mass moment of inertia of blade 
relative to drag hinge 
mass moment of inertia of fuse- 
lage about X and Z axis through 
A 

mass moment of inertia of shaft 
about X and Z axis through A 
coupling spring coefficient be- 
tween fuselage and shaft 
distance from the axis of rota- 
tion center to drag hinge center 
hydraulic damping moment coef- 
ficient 

nondimensional hydraulic damping 
coefficients 

dry friction moment coefficient 

nondimensional dry friction 

damping coefficients 

mass of the blade 

number of blades 

damping coefficient of blade 

radius of the rotor blade 

static mass moment relative to 

drag hinge 

time 

fuselage damping coefficient 
natural frequency of fuselage 
hinged mass ratio 
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damping coefficients 
natural frequencies 
hinged mass ratio 
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(C-model) 

fuselage damping coefficients 
natural frequencies of fuselage 
hinged mass ratio 
nS 2 H 2 /2I(J 0 +nm b H 2 ) 

angular motion coordinate of 
fuselage about X and Zi axis 
angular motion coordirate of 
shaft about X and Z a/.is 
angular amplitude of fuselage 
about X and Z axis 
angular amplitude of shaft 
about X and Z axis 
angular orientation of different 
rotor blades 

angular deflection of k-th blade 
natural frequency of rotating 
blade relative to drag hinge 
distance of center of gravity of 
the blade from drag hinge 
coordinate of the center of gra- 
vity of the blade system in 
fixed reference frame 
n ondimensi onal blade parameter 

J L v. h S/l 

rotor speed 


Introduction 

Helicopter ground resonance is a 
self-excited vibration phenomenon rather 
v »an a forced vibration. The safety 
standards for helicopter strength require 
that a helicopter with an articulated 
rotor must have enough stable margin when 
it is in contact with the ground, 
especially, in case of rough landing. Such 
standards are based on linear theory, but 
in fact the most common design of the 
landing gear has nonlinear damping 
characteristics. How much of an effect on 
ground resonance will be induced? Does the 
nonlinear damping characteristic tend to 
increase stability or not? These questions 
are waiting for investigators to solve. 
The present authors [1,2] and Tongue [3] 
have engaged in such research. The purpose 
of this paper is to further study such 
problems. 

The nonlinear damping in the landing gear 
of a helicopter brings into action both 
roll and pitch motion of the fuselage. In 
a series of publications [3-13] about 
helicopter ground resonance instability, 
this aspect was not taken into account by 
many authors. Also in the past, many 
investigators have used Coleman* classical 
multi-blade coordinate theory to simplify 


n x , nz 
p x» p z 

t c 

e, ,e, 

fx fz 
®sx' 0 sz 
0 fxO' e fzO 
0 sxO' 0 szO 

Y\ 

?k 

<*>s 

Pc 

V' $ 


'0 

a 


*»5 


PRECEDING PATS BLANK NOT FUMED 


the study of the rotor blade dynamics. 
Most previous authors have considered 
linear dynamic models. 


In the present paper, an approximate 
method has been used to calaulate the 
limit cycle behavior of a given helicopter 
model which has a nonlinearity in the 
landing gear such as a hydraulic 

resistance which is approximately 
proportional to the square of velocity or 
dry friction resistance. An experimental 
investigation has been carried out as 
well. The test apparatus is an improved 
version of Bielawa’s rotor model [4]. The 
experimental results are in good agreement 
with those of the theoretical analysis. 

In the present paper, we study the effect 
of combined roll and pitch fuselage 
motions, and consider the relationship 

between the lagging motion of the blades 
and the fuselage motion. Also discussed 
are the differences in the unstable 
boundary between the linear and the 
nonlinear model. We believe that these 

results will be helpful in further 
understanding the physical essence of 
helicopter ground resonance instability, 
and in protecting against such an 

instability. 


A pproximate method of nonlinec analysis 


The analysis is based on a helicopter 
model with an articulated rotor. A sketch 
of this model is shown in Fig.l. We will 
assume that the blades in the plane of 
rotation are rigid but articulated; the 
spring and damping coefficients of the 
fuselage and the shaft are isotropic, but 
their mass moments of inertia are 
anisotropic. Nonlinear damping includes 
hydraulic resistance which is proportional 
to velocity squared and also dry friction 
resistance; these can.be represented by 
M h ieie and M ,3 6ign0, respectively. The 
equations of motion for the model are six 
nonlinear differential equations with 
unknown functions 7 , « , 0f x , e fz * ®sx * nd 
B sz . These are as follows (see Ref.[l]): 
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In order to simplify the problem, we will 
discuss two special cases: 


1. A model with a single nonlinear 
damper. 

We suppose that the motion in the X 
direction is constrained, i.e, 0 SX «8 fx =0, 
and the subscript z in 0 sz and 0 f z is thus 

eliminated, and define that: 

§ s =m hz 9 s * ®f =m hz e f 
J H m hz SH / J sz' «=^hz SH/J xz 
Tiiut (1) becomes: 
ri +2n b >i +«> 2 >f-2a(?+n b ? ) =£ b § s 
? +2n b ? +u) 2 J+2Q(^+n b ^) =0 

y 2n 2 v< p ^v 2 n 4 V p Jv? (2) 

3f + 2 n69f+ (p2 + p2) §f + ^^sign§f 

- 2n 8§ -p 2 9 3 =0 

where *=m hz m d2 

For the nonlinear terms Of I Of I and signOf , 
we apply an approximate procedure which is 
called quasi-linearization. 

Let: Of = i Of I cos*t 

Representing the above nonlinear terms by 
a Fourier series and retaining the first 
harmonic produces (see Ref . [3] , [ 14] ) 
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Substitute these expressions into Eqs.(2), 
and require the determinant of 
coefficients to be zero for nontrivial 
solutions. The final characteristic 
equations are: 

real part: 
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imaginary part: 
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The coefficients a,b,c...etc. of the above 
Eqs. are given in Appendix 3 of Ref.[l]. 

For the case of Ksf-oo , <L-0f= 9. i.e the 
shaft is connected together with fuselage, 
the differential equations become a 
simpler nonlinear mathematical model for 
analysing helicopter ground resonance. 

"n + 2 +W g i[-2 Q( ?+n b 1 ) « £ a 5 

? +2n^-HoJs+2n(^+n b T)»o (5) 

0 + 2n Q 9 + 0 I 0 I + e = 

where : 

VW VW VW 2=Q/P 0 

For this case, the final characteristic 
equation is the same as in Ref [3]. 


Eqs. (6) is a set of nonlinear differetial 
equations for determining four unknown 
functions , C , 9 X and 8 Z . 

Using the method of harmonic balance and 
retaining only the fundamental harmonic, 
we take 

®x*®xO cosXt 

®z“®zc COSAt+ ® zs sinXt= e z 0 cos(xt-y 2 ) 

cosAt+7 g sinAt »? 0 cos(At-<^) t7> 
S*S C cosAt+f g sinAt =? 0 cos (At- <f f ) 

Here <f z , <f n , (f t aje respectively the phase 

difference of B , ? and ? with respect to 9 
z 1 x 

Mow we consider the term signfS^. and I ^ I . 
As was mentioned previously (2), we can 
obtain the following approximate result: 


Eq.(3]_ or (4) is eighth order in a. Also, 
the 0fQ angular amplitude nov; appears 
explicitly. This means that the limit 
cycle amplitude of the nonlinear system is 
directly related to A , the frequency of 
oscillation, and both must be solved for 
simul taneously . 

The Broyden method [15] was used to solve 
the above simultaneous equations. Let 
functions f ( Sf o » A ) an d f2(®fO» A ) he equal 
to the left hand sides of Eq.(3) and 
Eq. (4) , respectively. 


For a, 5 z g>0 
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Inserting Eq. (7) into Eq.(6) we obtain a 
system of eight simultaneous nonlinear 
algebraic equations for determining the 

quantities, n c .n s *?c« ?s> 8 xO'®zc< S zs A • 


fl<9 fo ,>0-0 

f 2 ( p , A ) =0 

2. A model with two nonlinear 
dampers. 

We suppose: 

K sf =0 ° ' e sx" e fx =0 x ' e sz =e fz =e z 
and define: 

J 0 =(J X +J Z )/2 ' *x =8 x' S z =e z' 1‘ISU/J o . 
?=?SH/J 0 . J x =J fx+ J gx , V J fz +J sz . 

W J x ' VV J z 


The general equations of motion for the 
model are given by Eqs. (6). 
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l^HxV^’xVx 1 ®x +m dx sign8 x = - £ 2 S 

V 2n *®z +P z 5 z +n h A 1 V m dz sign ®z* f 3 \ 


By elimination of the first four unknowns, 
the equations can be rewritten more 
compactly as a matrix equation with four 
unknowns, e x0 , § zc « ®zs* A • 

Let: 
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a= (u)2/^-i)/C l 
b=2n b /C 1 A 
c=2Qn fa /C ^a 2 
d=20/£ 1 A 

In order to solve 
IP) must be zero. 
Broyden's method 
equation. 

Exper i mental 



Eq . (8) , the elements of 
Here we also apply 
to solve the above 


Investigation 


The stability test 
disturbance, so two 
mounted on the test 
used to generate a 
direction of roll 
fuselage. 

The device 
calibration 
motor 
provides 
cam which 
that 
voltage 
actual 


needs an initial 
electro-magnets were 
equipment. Those are 
disturbance in the 
or pitch of the 


for angular amplitude 
is shown in Fig. 6. An electric 
with stepless variable speed 
a vibration source, and drives a 
can generate a sine wave, 3 uch 
the relationship between output 
from the transducers and the 
angle can be determined. 


1. Experimental equipment 


3. Frequency response test 


An overall view of the test model and 
recording equipment used in this study is 
shown in Figs. 2 and 3. The details of the 
various components of the model appear in 
Ref.[l] and [2]. The A-model has only one 
degree of freedom for fuselage motion. The 
B-model has two degrees of freedom for 
shaft and fuselage motion. The C-model has 
two degrees of freedom for the roll and 
the pitch motion of the fuselage. The 
characteristic curves of the hydraulic and 
dry friction damping are shown in Fig. 4 (a) 
and (b). The method of determining these 
parameters is described in Ref.flJ. 

2. Measurement system 

The X and Z< angular displacements and 
velocities of the fuselage and the shaft 
are obtained by WDT, R3 0D, velocity 
transducers located near the dry friction 
dampers, and a LVDT located near the 
gimbal support assembly. The output 
voltage of the transformers and velocity 
transducers is proportional to angular 
displacement and velocity of the roll and 
the pitch of the fuselage and the shaft. 
The output from the transformers are 
amplified and recorded on a multiple 
channel tape recorder, HP3968A. 

In order to obtain a one per revolution 
signal and the lagging motion of the 
individual blades during rotation of the 

blades, we mount a 13-ch brush and slip 
ring assembly on the shaft between the 
electric motor and hub, see Fig. 5. One of 
these slip rings is not a closed ring. It 
gives an impulsive signal once per 

revolution, so a very accurate measurement 
of rotor speed can be provided. Three 
angular transducers, R30A, are mounted on 
the drag hinges. The signals for lagging 
motion of each blade were amplified and 
recorded through brush and slip ring 
assemblies. The 8 channel signals were 
recorded simultaneously on a tape 

recorder, and further analysed by a 

Frequency Spectral Analyzer, HP3582A. 

Finally, the phase plane plots and 

frequency spectral plots were plotted by a 
X-Y recorder. 


The purpose of the frequency response teBt 
is to determine the natural frequencies of 
the system, and to provide an independent 
check on the system parameters. Before the 
test, the drag hinges are fixed, and a 
known block mass (3 oz) is put on a drag 
hinge. A centrifugal force caused by the 
bias mass will excite the system, and the 
frequency response versus rotor sneed can 

be found. 

4. Stability test 

Testing the model for self-excited 
instability regions was accomplished by 
slowly varying the rotor speed until 

instability was observed in response to an 
initial disturbance. 

Results and Di scussion 

Fig. 7 shows the limit cycle behavior when 
only the hydraulic nonlinearity is present 
in the A-model. The figure shows that the 
maximum limit cycle amplitude occurs near 
the critical rotor speed of the linear 
system, and that the dominant response 
occurs in the region, (1.7-2.0)*2 . 
Furthermore, there is an abrupt change in 
response from ( 1 . 9-1 . 95 ) *2 in Fig. 7. This 
means that the amplitude response is very 
sensitive to small changes in rotor speed 
in this range of rotor speed. 

Fig. 8 shows the limit cycle behavior with 
both hydraulic and dry friction 
nonlinearities. The curve has a shape 
similar to that of Fig. 7, but the response 
amplitude is smaller than that in Fig. 7. 
Physically the action of dry friction is 
to increase the equivalent viscous damping 
in the landing gear. The instability 
region becomss narrower, and the response 
amplitude reduces. 

Fig. 9 shows the limit cycle behavior of 
the B-model. There are two instability 
regions in the figure. The first one is 
dominated by the shaft motion. Because the 
effect of nonlinear damping of the 
fuselage on the Bhaft is weak, the shaft 
has a large limit cycle amplitude (for 


H8 


some range of a beyond the measurable 
angle) . In contrast with the first region, 
in the second instability region there is 
not a large limit cycle amplitude because 
of the direct effects of nonlinear damping 
of the fuselage on the fuselage motion. 

As far as the C-model is concerned. Fig. 10 
shows the results of frequency spectral 
analysis of the fuselage motion at a rotor 
speed of 1.2(HZ). The mode shape of the 
unstable motion is dominated by roll of 
the fuselage, and the oscillation 
frequency is equal to the roll natural 
frequency of the fuselage. Simultaneously, 
there is a small limit cycle amplitude in 
the pitch motion; its oscillation 
frequency is not equal to the pitch 
natural frequency, but is equal to the 
roll natural frequency. The phase angle 
between roll and pitch is equal to 51“ 
corresponding to a frequency of 0.84 (HZ). 
The phase curve in the vicinity of 
0.84 (HZ) is almost parallel to the 
frequency axis, and the phase angle is 
quite insensitive to frequency. 


conclusion which is derived by classical 
linear theory. However, the individual 
blades have different amplitudes and phase 
angles due to the actual small different 
blade damping characteristics. 

The equation of motion of the k-th blade 
can be written as follows: 

^k +2n ^k +W s r k = I^ ( * si"VZcosV 

The periodic term in the above equation 
can be removed by using the Coleman 
co-ordi nates. 

T 3 ZI fy. sin K 

k=l K * 

k-1 K k 

If it is expected that 1 and ? are well 
represented by a fundamental harmonic with 
oscillation frequency, X , then 1. may be 
defined as: K 


Fig. 11-12 show theoretical and 

experimental results of the limit cycle 
amplitude and phase angle behavior. A 
reasonable experimental verification of 
the theory is achieved. Because of the 
experimental complexity and calculation 
approximation, the quantitative accuracy 
of the verification is not very high. 
However, the general features of the 

experimental results are well predicted by 
the theory. For reference, in Fig. 13 the 
stability boundary results from linear 
theory are shown and compared to 

experimental data obtained with the 

nonlinear dampers removed. The agreement 
between theory and experiment in this case 
is good. 

Fig. 14 shows a phase plane plot which is 
directly plotted by a X-Y recorder. The 
measurement point of displacement and 

velocity is located at the same position. 
Their phase difference approximates to 90 s , 
the pattern displays an ellipse 
approximately, and the higher order 
harmonic components are included in the 
ellipse curve. From this pattern, we also 

can see the transient process leading to 
the steady state limit cycle amplitude. 
Different disturbances have different 
unstable processes. In this case the roll 
oscillation of the fuselage tends to a 
limit after about 7 pseudo cycles, and the 
transient process includes more higher 
order harmonic components. 

Fig. 15 shows the results of a frequency 
spectral analysis of the lagging motion of 
the individual blades at a rotor speed of 
2.04(HZ). It is found that in the 
instability region the individual blades 
have an identical oscillation frequency, 
p, which is equal to the difference 
between the rotor speed and the 
oscillation frequency, a. , of the fuselage. 
This result is in agreement with 


-^k = V os < Pt+27 ’'*/n ) 

i.e the following conditions must be 
satisfied: 


* P= Q - 

* the amplitudes of the individual 
blades are identical. 

* the phase angle of the k-th blade is 
eq ual to 2 k */n . 


In the present experiment, only the first 
condition is satisfied. The modest 
violations of the latter two conditions 
are the main reasons for the error between 
theoretical and experimental results. 


In order to provide a better physical 
understanding of the unstable motion, we 
discuss the relationship between motion of 
the hub center and the center of gravity 
(C.G) of the blades. In the fixed 
coordinate system XOZ, the motion of the 
common center of gravity of the blade 
system is expressed by 


-v x ~ £ h sin \ 

<=! * 

Z =z+ — y^f.cosi/: 
n k=l k 

X= h£_ 8 cosA.t 

i=l xOi i 


i-1 


8 n . cos (A. . ) 
zOi i J zx 


Fig.lfi shows the time history of and S . 
The test curves include components with 
frequency (22-x), and a , i n addition, a 
cos xt component. This also can be clearly 
found from Fig. 17. The motion locus of the 
C.G of the blades is not exactly an 
ellipse, but a somewhat more complicated 
curve. The starting point at t-0 does not 


coincide with the end point at t»2(sec.), 
howevec over a much longer time interval, 
the motion does appear periodic. The phase 
angle between hub center and C. G of blades 
is equal to 46* from calculation, and 15? 
from test. The phase angle depends upon 
the dampings of the fuselage and blades, 
as well as rotor speed. 

Fig. 18 shows the theoretical and 
experimental results for a rotor speed of 
1.2(HZ). The test curves are approximately 
cosine waves. The motion locus of the C. G 
of the blades approximates to an ellipse, 
and excellent agreement of calculation 
with test is achieved for the motion locus 
of the hub center. As compared with Fig. 17 
and 18, we find that the difference 
between theory and experiment for the hub 
center motion mainly comes from the 
variation of blade motion from one blade 
to the next. 

Fig. 19 shows the test results for a rotor 
speed of 1.96 (HZ). All higher order 
harmonic components are included in the 
motion of the fuselage and blades. Their 
motion loci have guite complex shapes, but 
they always tend to a limit amplitude. 


Coaclufliona 

A number of conclusions may be reached. 

1. Unstable motion of a helicopter in 
ground resonance can occur only when the 
rotor speed is within a certain region, 
and at the same time a initial disturbance 
acts on it. This conclusion is suitable 
for both linear and nonlinear systems. 

2. During unstable motion the 
oscillation frequency is equal to the 
natural frequency of the linear system, or 
approximately the natural frequency for a 
nonlinear system. 

3. The nonlinear dynamic rotor system 

with multiple nonlinear dampers has more 
complex motions than the system with a 
single nonlinear damper. The coupled 
notion of roll and pitch of the fuselage 
is due to blade lagging motion. This 
coupling motion is important in 
understanding ground resonance of an 

actual helicopter. 

4 . The agreement of theory with the 
experimental results verifies that the 
analytical method which is described in 
this paper is useful and reliable. 

5. The different characteristics of 

individual lag blade dampers will lead to 
different limit cycle behavior. 

6. The results show that a landing 

gear with hydraulic nonlinear damping 
characteristics and dry friction damping 
may be advantageous for protecting against 
ground resonance instability. It is 


unlikely, perhaps, that a helicopter 
resigner would rely on a nonlinear damper 
for primary protection against ground 
resonance. However, such a damper may be 
useful as a failsafe device to ensure a 
benign failure mode should a primary 
linear damper become inoperable or prove 
inadequate under some operating 
conditions. 

7. It is expected that the design and 
operation of blade dampers will also be 
very important for avoiding ground 
resonance. This is a logical topic for 
future study. 
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Jerry Mlao, Sikorsky Aircraft : I would like to make a comment. Earlier Jing Yen said tha* - we 
may not be able to analyze ground resonance, but apparently from this paper we can see with a 
nonlinear damper we can even calculate the amplitude of the oscillation. Therefore, if you orJy 
know what numbers to put into the analysis you will be able to calculate it. Thank you. 
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